Nuclear Spin Relaxation Rate of Disordered px + ipy-wave Superconductors 
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Based on an effective Hamiltonian with the binary aUoy disorder model defined in the trian- 
gular lattice, the impurity scattering effects on the density of states and especially on the spin- 
lattice relaxation rate l/Ti of Px + ip^-wave superconductors are studied by solving numerically 
the Bogoliubov-de Gennes equations. In the clean limit, the coherence peak of 1/Ti is observed as 
expected. More intriguingly, for strong scattering potential , the temperature dependence of 1/Ti 
exhibits the two different power law behaviors near Tc and at low temperatures, respectively, which 
is in good agreement with the nuclear quadrupolar resonance measurement. 

PACS numbers: 74.20.Rp, 74.25. Jb, 74.25.Nf 



The novel superconductor, Na2:Co02:yH20 [x = 0.35), 
recently discovered by Takada et al [l| has stimulated 
much theoretical and experimental interest in studying 
its spin and orbital symmetries of the Cooper pairs, which 
helps to explore and understand the underlying supercon- 
ducting mechanism of this new material. At present this 
issue is hotly debated and there still exist controversies 
both theoretically and experimentally. Resonating va- 
lence bond 0] theories of the triangular lattice t-J model 

Q 111 describe Na2,Co02-yH20 as an electron-doped 
Mott insulator based on the fact that the Co**"*" atoms 
has spin-i as Cu^+ in high-Tc cuprates. Such theories 
prefer the spin-singlet c? + id'-wave pairing {d = dx2^y2 
and d' = dxy)- On the other hand, theories based on 
a combined symmetry analysis with fermiology 0] favor 
triplet over singlet pairing with the d-vector perpendic- 
ular to the cobalt plane [most probably d — z{px ± ipy), 
i.e. Px ± ipy-waye pairing state]. As for the experi- 
mental aspect, a series of nuclear magnetic/quadruplolar 
resonance (NMR/NQR) measurements i J., also 
give divergent implication about the pairing symmetry in 
Naa;Co02-2/H20. Namely, for spin pairing Ref. ^ found 
the signature for the triplet pairing while Ref. ^ sup- 
ported the singlet pairing although both examined the 
temperature dependence of the Knight shift; for the or- 
bital wave function of the Cooper pairs, a fully-gapped 
superconducting state 0; Hi is inferred from the exis- 
tence of the Hebel-Slichter coherence peak of the 
spin-lattice relaxation rate (1/Ti); on the other hand, 
opposite conclusion is drawn from the NQR experiment 
[lO] that l/Ti decreases across without the coherence 
peak and follows a power-law dependence deviating the 
exponential relation, which implies a nodal gap function. 

Motivated by the novel power law behavior reported in 
Ref. 01 and the fact that various impurities or defects, 
such as inhomogeneity in the Na content, oxygen vacan- 
cies and substitution atoms in the C0O2 layer etc., are 
present during the synthesis and handling of this com- 
pound, in this Letter, we investigate and elucidate the 



effect of disorder on the nuclear spin relaxation rate of 
the Px ± ipy-wave superconductor for the first time. Our 
results indicate that the seemingly incompatible exper- 
imental observations of the temperature dependence of 
l/Ti may be reconciled within the picture of the disor- 
dered chiral px + ipy-wave pairing state if the effect of 
impurities is properly considered. 

Here we employ a mean-field Bogoliubov-de Gennes 
(BdG) Hamiltonian on a tight-binding triangular lattice 
with the nearest-neighbor (NN) hopping integral t and 
the NN bond pairing potential Aij resulting from the ef- 
fective attractive interaction V. The model Hamiltonian 
is expressed as 



Hp 



J2[A,,{cl^cl+clc]^)+h.c], (1) 



where tij = t is the NN hopping integral and in the re- 
minder of this paper, we choose t < according to the 
analysis on the band calculation |0| and the energies 
will be measured in unit of jtj. /i is the chemical po- 
tential. Here, we adopt the binary alloy disorder model 
[1^ where ei is the (5-function-like scalar scattering po- 
tential and takes the value Uq on certain lattice sites ac- 
cording to the impurity concentration ni,„p and zero else- 
where. The spin-triplet pairing potential Ay is defined 
as Aij = Y{{cii;Cji) + (cijCj-f)). Note that only the dz 
component of the spin-triplet pairing is considered here 
in view of the experimental indication that the d vector 
is parallel to z-axis according to the invariant behavior 
of the Knight shift for the in-planc magnetic field jSj • bi 
the homogenous case, the px + ipy-wave pairing state is 
expressed as 

Ap^+ip^{k) = 2Ap sm{kx) +sm{kx/2)cos{Viky/2)+ 
iV3cos{kx/2)sin{^/3ky/2)\ (2) 
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where Ap = J2i,5 Au+se'''^'^^^ with i + 5 the six NN 
sites of i. By applying the self-consistent mean- field 
approximation and performing the Bogoliubov transfor- 
mation, diagonalization of the Hamiltonian H^ff can be 
achieved by solving the following BdG equations: 



j 



Er, 



(3) 



where u", v" are the Bogoliubov quasiparticle amplitudes 
with corresponding eigenvalue Hij = —tij -f ^ij (ei — 
^). is calculated according to: A^j = -j '^n^u^^^j* + 
tanh(2^y)- Throughout this work, we set V = 
2.3 and fi = 1.0 which gives rise to Ap = 0.12 and the 
electron number per site is 1.39 in the absence of disorder. 
Due to the vanishingly small anisotropic factor of the 
Px + ipy-wave pairing 115|, there is one s- wave-like full 
gap opened at approximately Acap ~ 0.4 (see Fig. ^ for 
the clean limit). 

Once the self-consistent quasiparticle spectrum is ob- 
tained, the nuclear spin relaxation rate is calculated ac- 
cording to [13 1 



n J j[pUE)pl{-E')-f^A{E)pl{-E')] 
xf{E)[l^f{E')]6{E-E')dEdE', (4) 
-ttT / p^,{E)p^^,{-E)f{E)dE (5) 
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where Pap{E) is expressed as 



p\\iE) p[^,{E) \ 
pl\{E) pgiE) 
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For the unconventional pairing, the contributions from 
the off-diagonal elements pi2 and p2i are absent in Eq. (O 
due to the pairing symmetry which forbids the on-site 
pairing amphtude 18]. The local density of states (DOS) 
is calculated according to 



p{t,E)^2pt,{E)^2p^,{~E). 



(6) 



To numerically investigate the disorder effect on the 
electronic structure and accordingly the nuclear spin re- 
laxation rate in the 2D system, for certain impurity con- 
tent nimp, typical DOS and l/Ti are obtained by averag- 
ing over 20 impurity configurations with the size of the 
supercell 20 x 20 and 100 wave vectors in the supercell 
Brillouin zone 14]. For each impurity configuration, self- 
consistent bond pairing potential A^ is obtained with 
the maximum relative error between two consecutive it- 
eration steps is less than 10"'^. In this work, assum- 
ing that 1/Ti(i) is determined by the on-site R{i, i) with 
ignoring the minor contributions from neighboring sites 
R{i,i + S), we have the impurity-averaged 1/Ti = Tj~^(i) 
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FIG. 1: Density of states as a function of energy of the dis- 
ordered Px + ipy-wave pairing state for the clean limit (solid 
line), riinip = 2%, Uq = 2 (dashed line), riimp = 5%, Uo = 2 
(dotted line), nimp = 5%, Uo = 10 (dash-dotted line). Also 
shown in the inset is the temperature dependence of the 
impurity-averaged order parameter Ap(r). 



and p{E) = p{i,E), where (. . .) denotes averaging over 
space and impurity configurations. 

In Fig. n] we illustrate the dependence of the disorder- 
averaged DOS on the impurity content n.imp and the scat- 
tering strength Uq. The set of curves display clearly (1) 
the shrinking (and even vanishing) of the energy gap (2) 
the smearing and decreasing of the coherence peak as 
nimp and Uq increase. As nimp increases from to 2% 
and 5% for the fixed Uq = 2, the gap is filled from the 
gap edge, resulting in a reduced effective gap. And when 
C/o — 10 for riimp = 5%, strictly speaking the energy 
gap is closed with finite p{E) a& E Q (although the 
DOS hump around Acap can still be identified) and the 
residual DOS at the Fermi level is as large as 60% of the 
normal state value. The large residual DOS in the su- 
perconducting state is consistent with the experimental 
study of specific heat indicating the importance of 
inhomogeneity in this material. In the inset of Fig. ^ 
we give Ap(T) in the clean limit, riimp — 2%, Uq = 2, 
nimp = 5%, Uq = 2 and rUmp = 5%, Uq = 10, showing 
that both the order parameter and the transition tem- 
perature are reduced and the decrease of Ap as T — s- Tc 
is also not so sharp in the disordered cases as in the dilute 
limit. These behaviors are in consistence with the pair 
breaking effect of non-magnetic impurities in unconven- 
tional superconductors. 

The evolution of the density of states originates from 
the impurities, and depends on their content as well as 
strength. For the weak scalar scatterer with Uq = 2 (cor- 
responds to Cs ~ 0.5 according to Ref. T^, there are two 
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peaks of the local DOS around an isolated impurity near 
the gap edge with energies Wb/Acap — ±0.89 highlight- 
ing the presence of impurity bound state in the uncon- 
ventional Pa;+ipy-wave superconductor 18, in remark- 
able contrast to the conventional s-wave superconductors 
although both are fully gapped. And for strong poten- 
tial such as Uq = 10 (corresponds to Cg — 2.5), we have 
Wb/Acap — ±0.37. For finite impurity density, to look 
into the closure of the energy gap induced by the impu- 
rities, we study the impurity-averaged self-energy of the 
quasiparticle, which is determined by the self-consistent 
equations 20, 21]: 

G{k, iujn) = [iwn - i{k)az ~ Ap^a-^ - ^PyC^y] (7) 
g{iujn) 



where iuJn 



Ump 2 If \ 



(8) 



and Ap^ are real 



and imaginary parts of Eq. 121, respectively, ^(fc) is 
the normal state quasiparticle energy and giiujn) = 
J (Pk/{2TTyG{k,iujn)- For the weak scatterers Uq^ '3> 
g{iujn), the scattering rate 7 [determined by — )■ 0) = 

-i^] is 7 = y {-nNom^pUlY - A^g^p, where Nq is the 
normal density of states per spin at the Fermi level. 
Therefore, the impurity parameter nimpt/o must be larger 
than Acap (r)/7riVo to entirely close the gap, i.e. 7 ^ is 
real. This effect of disorders for the Px + ipy-wave pairing 
is significantly different from that for the high-Tc nodal 
d2.2_j,2-wave pairing. In the d-wave pairing, infinitesimal 
'T-impC^o gives rise to finite density of states at the Fermi 
level. In the strong scattering limit, Uq^ <C g{iujn), we 
obtain 7 = niinpAcap/TriVo, which is the same as the 
result of d-wave superconductors with unitary impurities 
[2^ . The above discussion qualitatively explain what 
we illustrate in Fig. ^ 

The impurity effect on DOS is manifested by the vari- 
ation of NMR relaxation according to Eq. 0. First, 
we address the temperature dependence of l/Ti in the 
absence of disorder (clean limit). When njmp — 0, the 

— Aq^^p divergence of the DOS p{E) at the gap 
edge Agap will lead to the Hebel-Slichter coherence peak 
of 1/Ti just below Tc as shown in Fig. [3 although the 
jump of the peak is much lower than that of the s- 
wave pairing because the coherence factor changes from 
1 -I- Agg^p/E'^ in the s-wave case to 1 in the Px + ipy-wave 
case. Furthermore, when T < Tc and Acap > T, l/Ti 
of both the s-wave and px + ipy-wave pairing states shows 
the behavior of e~^'^''p/"^ due to their fully gapped na- 
ture. Also shown in Fig. |21is the behavior of 1/Ti of the 
gapless px-yfave pairing state for comparison. As for the 
nodal p^-wave pairing [real part of Eq. Q], the slower 
logarithmic divergent of p{E) and the halved coherent 
factor result in much suppressed coherence peak and in 
the low temperature region p{E) cx E gives rise to the 
dependence of the nuclear spin relaxation rate as shown 
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FIG. 2: Temperature dependence of l/Ti normalized by 
its value at Tc for the px + ipy-wave (solid line), s-wave 
(dashed line), and p^^-wave (dotted line) superconductors. 
Solid squares represent the experimental data extracted from 
Fig. 2 of Ref. Efl 
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FIG. 3: Temperature dependence of Ti/Ti for the disordered 
Px + ipy-wave when nimp ~ 5%, Uo = 10 (solid line) and 
?T-imp = 5%, Uo — 2 (dashed line). The Solid squares are the 
same as those of Fig. H 



in the figure. 

Figure |31 shows the nuclear spin relaxation rate when 
the impurity scattering is present. Tw o typ ical results are 
compared with the experimental data[l£j. As expected, 
the coherence peak of 1/Ti disappears for both the weak 
and strong scattering cases, in accord with the finding 
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that small concentration of impurities is able to smear 
the sharp divergence of the DOS near the gap edge. For 
strong disorders with Uq = 10, our results of 1/Ti as a 
function of temperature indicate a T^'^ dependence at 
the low temperature region and a T^'^ relation near Tc, 
being consistent with the results of an approximate con- 
stant p{E) near zero energy as shown in Fig. ^ More 
importantly, these results are in good agreement with 
the experimental observation of Ref. On the other 
hand, the study of weak disorders with Uq = 2 exhibits 
that 1/Ti first drops with a T" (n ~ 3) law at the vicin- 
ity of Tc and then exponentially similar to its behavior in 
the clean limit showing a downward curvature. The oc- 
currence of the exponential behavior is due to the open- 
ing of the energy gap at certain temperature T* for the 
weak scatterers. And T* is governed by the solution of 
AGap(T*) = 7riVo"imp?7o' which results in T*/Tc ~ 0.8 
for Uq = 2 and riimp = 5%. According to this, one will 
have a power-law dependence of 1/Ti down to lower tem- 
perature region by simply increasing the impurity content 

^imp • 

In summary, we have for the first time elucidated 
the disorder effects on the electronic structure and nu- 
clear spin relaxation rate of the Px + ipy-wave pairing 
state, which is closely relevant to the new supercon- 
ductor Nao.35Co02-yH20. The experimentally observed 
temperature dependence of l/Ti is explained satisfacto- 
rily. It is also interesting to compare the present re- 
sults with those for another layered superconductor , 
Sr2Ru04|3 with a possible px + ipy-wa.ve pairing 
symmetry. We found that the gap in Nao.35Co02-yH20 
is highly isotropic |l5| in the triangular lattice, while the 
gap in Sr2Ru04 with the proposed Px + ipy-wa.ve pairing 
is strongly anisotropic l24j. Therefore, we predict that 
future l/Ti measurements on high-quality samples of 
Nao.35Co02-j/Il20 may observe both the coherence peak 
and the exponential dependence on temperature, in con- 
trast to the T'^ behavior of l/Ti in the significantly 
anisotropic Sr2Ru04 even if both are in the hypothetical 
Px + ipy-'wsLve pairing states. 

On the other hand, as far as the DOS-related physical 
quantities, such as 1/Ti, are concerned, there should be 
no qualitative differences between the Px + ipy-wave pair- 
ing and the d + id'-wave pairing state in the sense that 
they are both fully-gapped and have similar response to 
impurities [l^. Therefore, more definitive experimental 
measurements of the Knight shift are demanded to de- 
termine the symmetry of the spin part of the Cooper 
pair wave function. Moreover, the muon-spin-relaxation 
measurement j26j , which is sensitive to the time-reversal- 
symmetry breaking effect in the px + ipy- and d + id'- 
wave pairing states, and the phase-sensitive Josephson- 
tunnelling related experiments can give more decisive 
evidences to distinguish the gapped pairing states from 
the nodal ones, such as the Px and c?j,2_„2 waves. The 



scanning tunnelling microscopy experiments |28l | are also 
able to shed light on the pairing symmetry by examining 
the energy and spatial pattern of the impurity and 
vortex 0,113 states. 
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